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Abstract. Second order linear non-autonomous diflfcrential equations with 
negative stiffness are considered. Using Chetaev-like (Lyapunov-like) func- 
tions, necessary (sufficient) conditions are found for the solutions to be bounded 
for all initial conditions if any one of the coefficients is constant. Conclusions 
are then extended to include systems where both coefficients are time- varying. 



1. Introduction 

The second order linear differential non-autonomous homogeneous equation 

X + ait)x - k{t)x = , (1.1) 

and its non-homogeneous version are usually found in mechanical problems like 
the study of vibrations, and in perturbative analysis of solutions to differential 
equations in Physics. Here, a, fc : I — *■ M, dot stands for differentiation with respect 
to the independent variable t E 1= [to,(X)). Following Mechanics notation, a{t) is 
going to be called the damping of the system, while —k{t) will be referred as the 
stijfness. 

Too often a k{t) strictly positive for all t > to is automatically and tacitly 
regarded as a signal of instability of the trivial solution of the above equation. The 
reason for that is just that the corresponding to Ea. ljl.ll) planar system, 



1 



k{t) - a{t) 



(1.2) 



(x"^ = 1x1X21, Xi = X, X2 = x), has eigenvalues X±{t) oc ~a{t) ± y/ a{t)^ + 4fc(t), 
pointing to a saddle-like instability for any given value of t. Thus, if these eigenval- 
ues are slowly changing the system is assumed to have no bounded (finite) solutions. 

However, a rigorous foundation of this folk theorem is still missing as well as 
a precise definition of slowly changing. A number of authors has addressed this 
problem (see for instance El El E]) but either they deal only with positive 
stiffness or their results are not expressed as direct conditions on the time-varying 
coefficients. 

The aim of this paper is to help to fill this gap in the analysis of the stability 
on non-autonomous linear systems. With that in mind, in the next section, the 
Lyapunov's direct method is used to state two theorems for the necessary and suffi- 
cient conditions on the time variation of Q;(t) and k(t), for Ea. (|l.l(l to have bounded 
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solutions for any initial condition xq = x{to). Each theorem corresponds to a case 
where one of the coefficients is time- varying while the other one remains constant. 
In the last section, corollaries and examples are presented allowing to extend the 
conclusions to cases where both coefficients are time-varying. Throughout this 
manuscript, it is assumed that conditions are imposed on H1.I|I such that existence 
and uniqueness of its solutions are warranted for all initial conditions xq G and 
for all initial times to. The Euclidean norm, || • ||, in is used. The definitions 
and theorems for boundedness and stability in the sense of Lyapunov can be found 
in many textbooks (see for instances 



2. Main results 

Theorem 2.1. Consider the linear non-autonomous equation 

x + f3x-k{t)x = 0, (2.1) 

where (3 is a real constant and k{t) is a differentiahle real-valued function strictly 
positive in I. Then, 

(!) /3>0, 
(2) kit) < 

are necessary conditions for the solution x = o/ ^2.1\) to be stable in the sense of 
Lyapunov. Moreover, 



(3) 



> - 



_k{t)_ 

along with the conditions above, is sufficient for all solutions to be bounded ^xq. 

Proof. Condition follows straightforwardly from requiring no continuous in- 
crease of the volume of the phase flow due to the linear vector field v = A{t)x. In 
other words T', divv = Tr^(t) < 0, where div stands for the vector field divergence 
and Tr for the matrix trace. If /3 = the solutions will be unbounded since even in 
the limit k{t) — s- +0 the solution blows up linearly. 

To prove condition J^J recall first that for k{t) = the eigenvalues of the corre- 
sponding planar system are A± oc — /? ± \/ (? ^- 4fc, and x = is indeed a saddle 
point. So, assume now k(t) > 0. Then, consider the Chetaev's function [1], 

V^{x,t)^^k{t)xl^^xl. (2.2) 

Let Di be the circle defined by ||x|| — e^, let D2 be the union of the two subspaces 
distributed along X2 — with borders 9I?2 given by X2 — iy^k(to)xi, and let 
D C Dl n D2/dD2. Then, for k{t) > 0, and e as small as desired, Vx G D: 
V^{x,t) > y^(x, to) > 0. V)^{x,t) — only for x £ dD2, what includes the origin. 
The time derivative of H2.2|l along the direction of the vector field v is 

^ = lk{t)xl+Pxl. (2.3) 

For any x & D and for all i e J: V^{0,t) = and V^{x,t) > 0. Therefore, any 
solution of (|2.1|) with /3 > 0, fc(<) > and initial condition Xq G D will eventually 
escape to infinity, independently of how close Xq is chosen to be to the origin. With 
these conditions, the trivial solution of H2.1|l is unstable. 
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The sufficient condition ^ can be proved using the Lyapunov's function 

VA{^,t)^^k{t)xl + ^xl. (2.4) 

For k{t) > this function is non- negative and is equal to zero only at the origin. 
Since k{t) is differentiable in [to, oo) and decreases with time, Vx 7^ and Vt G I : 
Va(x, io) > V^(x, i) > 0. Va(x, io) 00 only when ||x|| 00. The corresponding 
time derivative along the direction of the vector field v is, 

^ = hit)xl + 2k{t)xiX2 
which can also be written as, 

ifcW k{t) 
k{t) - (3 



dVA 
dt 



(2.5) 



(2.6) 



This function is negative definite if Vt > Iq: k{t) < 0, /? > and -k{t)/k'^{t) > 2/13. 
With these conditions, solutions of system (|2.1() will be bounded Vxg. □ 



Theorem 2.2. Consider the linear non- autonomous equation 

X + a{t)x — uP'x — , 



(2.7) 



where lo is a real constant and a{t) is a differentiable real-valued function in I. 
Then, 

(1) a{t) > 0, 

(2) a{t) > 0, 

(3) Alna{t)>^[^+u;f 

Vt > to are together sufficient conditions for solutions of \2. 7| ) to he hounded Vxq . 
Proof. This theorem can be proved using the Lyapunov function. 



^4(x,t) = - 



1 1 



2a{t) 



(2., 



and the corresponding time derivative along the direction of the vector field v, 



dVx T 
— ;— = X 
dt 



2 dt 







1 


to 


i"(t)J 




2 


a.{t) 








~a{t) 



(2.9) 



□ 



3. Corollaries and examples 

Corollary 3.1. There exist time-varying coefficients ol{t) and k{T) such that so- 
lutions to Eg. Hl.lp will he hounded ^^.q. 

Proof. 

Example 3.2. The non- autonomous equation 

x" + a(r)x' - fc(T)a; = , (3.1) 
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where prime stands for derivative with respect to t, and Vr > tq 

k{T) = (3.2) 
= /3/(r)t-|/(r)-V(r)', (3.3) 

for all real n > 0, /3 > 0, /(r) > 0, /(r)' < 0, and /(r)- 1 [l//(r)]' > 2wV/3. 
To test this, change the independent variable in Ea. H2.1|l using dr/dt — /^"/^. 

□ 

Corollary 3.3. In Eg. J^l.l]) the time-varying coefficient a does not need to he 
strictly positive for the solutions to be hounded Vxq. 

Proof. 

Example 3.4. The non- autonomous equation 

x" + a{T)x' - k{T)x = Q , (3.4) 

where Vr > tq, and for all real lj , n > 0, [3 > Q, 

k{T) = u:^9{Ty\ (3.5) 

a{T) = /35(T)^+^-f5(r)-\9(r)', (3.6) 

*/5(t) > 0, g{Ty > 0, and d\ng{T)/dT > ^(t)"/^ [/3-i5(r)-i + c.2]2/2. 

And this can be tested using the parametric change dr/dt = g^"/^ in Ea. H2.7|l . 

□ 

Intuitively, a positive damping can control the instability induced by a negative 
stiffness if the former is asymptotically larger than the absolute value of the last. 
Theorems 12.11 and 12.21 together with corollaries 13.11 and 13.31 give a notion of how 
large the ratio \a{t)\/k{t) must be at any given time for the corresponding solution 
to be bounded Vxq . 

Remark 3.5. The sufficient condition \2. 11!^ for the existence of hounded solutions 
of Eq. \2.1]) is rather weak. 

Indeed, if the stiffness decays with respect to some finite positive damping (3 the 
solution of Ea. H2.1|l will asymptotically look like Ci + C2 exp(— /3t/2). For instance, 
for equation. 

Example 3.6. 

with TO > 1, condition \2. ll'A reads 

TO p 

So, the solutions are hounded even for m — 1 if P > luP' . For to > 2 the condition 
is satisfied asymptotically for any (3 and lo. 

For m — 2 the corresponding solutions are given by 



^ /7 _/lt^/\/l + 4w2 i3t\ / Vl + 4cj2 f3t\ 
r{t)^CiVie =*I ,^\+C2Vie 2*/c ,^ , 
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where C'l and C2 are the integration constants, and X(a, z) and /C(a, z) are the mod- 
ified Bessel functions of the first and second kinds, respectively. In the asymptotic 
regime, this solution can be approximated by 



(3.8) 



converging uniformly with i ^ 00 to the constant given by the initial condition xq- 
For completeness, note that in the limit of strong decay of the negative stiffness. 



Example 3.7. 



(3.9) 



the solutions for m > are given by, 

x{t) ^ Cie-^tj ( 



m m 

which for very large t can be approximated by 



/3 2we--' 



m 



m 



1 



/3 



+ C2 



m 



-fit 



converging uniformly with < — > 00 to the constant given by xq. 

Remark 3.8. Conversely, the sufficient condition \2. UtA for the existence of bounded 
solutions of Eg. ^^7f\j is very strong, it is close to necessary. 

For instance, for equation, 

Example 3.9. 



X + pf^x - iJ'x = , 
with TO > 0, condition \2. reads 

-I 2 



(3.10) 



2 

TO 



/3 



>0, 



which is not satisfied for any real (3 and uj if t 00. 

In fact, for < m < 3, Ea. H3.1()|l can be solved in term of Heun series, each one 
diverging for all positive /3 and real uj. Nevertheless, this condition is not actually 
necessary since in the limit of strong growth for the damping, equation 



Example 3.10. 

has solution, 



X + /3e"*i - w^a; = , 



(3.11) 



\ z TO 2 m / \2 m 2m J 

+ C2e^H-^<="') 



UJ 1 (3 



UJ 1 (3 



^ o > -^e™' - /C - + -, --e' 

\2 m 2m J \2 m 2m 

for any real cj, /? > and to > 0. Asymtotically it behaves like 

C,-^- + C.V^^-e~^^ , (3.12) 

also converging uniformly with < — > cx) to the constant related to xq. 

Finally, note that substitution dr/dt — t (with tq > 0) will change Ea. H3.9|l into 
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Example 3.11. 

x" + {P+l)^x' -u^^x^O, (3.13) 

and Eg. H3. 10(1 respectively into 
Example 3.12. 

x" + ipT"' + l)-x' -uj^^x^O. (3.14) 

The last two examples are generalizations of the Cauchy-Euler equation. Dif- 
ferently from the Cauchy-Euler case (to = 0), these generalizations have bounded 
solutions Vxo, provided the conditions discussed in this paper are satisfied. 
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